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Abstract 

We generalize Wigner's causality bounds and Bethe's integral formula for the effective range to 
arbitrary dimension and arbitrary angular momentum. Moreover, we discuss the impact of these 
constraints on the separation of low- and high-momentum scales and universality in low-energy 
quantum scattering. 

PACS numbers: 21.45.-v, 21.45.Bc, 34.50.-s, 34.50.Cx, 03.65.Nk 



1 



In quantum mechanics causality requires that no scattered wave propagates before the 
incident wave first reaches the scatterer. This constraint ensures a scattering amphtude 
which is analytic in the upper half plane as a function of energy B]- For the case 

of finite-range interactions the constraints of causality were first investigated by Wigner 
The time delay between the incoming wave and the scattered outgoing wave can be 
computed from the energy derivative of the elastic phase shift, At = 2hd5/dE. If dS/dE 
is negative, this implies a time advance of the outgoing wave. However, the time advance 
cannot be arbitrarily large since the incoming wave must first enter the interaction region 
before the scattered wave can exit. Since the derivative of the phase shift with respect to 
the energy is involved, this argument places a bound on the effective range of the scattering 
amplitude. 

Wigner bounds are particularly interesting in the context of low-energy scattering and 
universality. Universality at low energies arises when there is a large separation between 
the short-distance scale of the interaction and the long-distance scales given by the average 
particle spacing and thermal wavelength. One example of low-energy universality is the 
unitarity limit, which refers to an idealized system where the range of the interaction is zero 
and the S'-wave scattering length is infinite. It has been studied most thoroughly for two- 
component fermions. In nuclear physics, cold dilute neutron matter is close to the unitarity 
limit. However, most recent interest in unitarity- limit physics is driven by experiments 
with cold ^Li and atoms using magnetically-tuned Feshbach resonances. For reviews 
of recent cold atom experiments, see Refs. 0, [sj. Theoretical overviews of ultracold Fermi 
gases and their numerical simulations are given in j^, 0]. A general review of universality 
at large scattering length can be found in [8!|. 

Several experiments have also investigated strongly-interacting P-wave Feshbach reso- 
nances in ^Li and ^''K [9-13|. An important issue here is whether the physics of these 



strongly-interacting P-wave systems is universal, and if so, what are the relevant low-energy 
parameters. A resolution of these issues would provide a connection between, for example, 
the atomic physics of P-wave Feshbach resonances and the nuclear physics of P-wave alpha- 
neutron interactions in halo nuclei. Some progress towards addressing these questions has 
been made utilizing low-energy models of P-wave atomic interactions and P-wave 

alpha-neutron interactions j20-23||. A renormalization group study showed that scattering 
should be weak in higher partial waves unless there is a fine tuning of multiple parameters 

In this letter, we answer the question of universality and the constraints of causality for 
arbitrary dimension d and arbitrary angular momentum L. Our analysis is applicable to 
any finite-range interaction that is energy independent, non-singular, and spin independent. 
We present generalizations of Bethe's integral formula for the effective range j25| and Wigner 
bounds for arbitrary d and L. Our results can be viewed as a generalization of the analysis 
of Phillips and Cohen [26[, who derived a Wigner bound for the S'-wave effective range for 
short-range interactions in three dimensions. Here we show that for 2L + d > 4, causal 
wave propagation produces a fundamental obstruction to reaching the scale-invariant limit 
for finite range interactions. Instead we find the emergence of the effective range as a second 
relevant low-energy parameter that cannot be tuned to zero without violating causality. 
For the case of shallow bound states, we show that this second low-energy parameter also 
parametrizes the size of the bound-state wavefunction. Complementary work was carried 
out by Ruiz Arriola and collaborators. A discussion of the Wigner bound in the context 
of chiral two-pion exchange can be found in 27|] while correlations between the scattering 
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length and effective range related to the Wigner bound were discussed in H . 

We consider two non-relativistic spinless particles in d dimensions with a rotationally- 
invariant two-body interaction. We let L label the absolute value of the top-level angular 
momentum quantum number [29, For d > 2, L can be any non- negative integer. For 
d = 1, the notion of rotational invariance reduces to parity invariance. Here we assume 
a parity-symmetric interaction and write L = for even parity and L = 1 for odd parity. 
We analyze the two-body system in the center-of-mass frame using units with h = 1 for 
convenience. With reduced mass and energy p^/(2/i), we rescale the radial wavefunction 



as 



{prf-'^"Rf,Ar). (1) 

The interaction is assumed to be energy independent and have a finite range R beyond which 
the particles are non-interacting. Writing the interaction as a real symmetric operator with 
kernel H^(r, r'), we have the radial Schrodinger equation 
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The normalization of ui'^yr] 



is chosen so that for r > R, 
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Here and Ya are the Bessel functions of the first and second kind, and Sl^p) is the 
phase shift for partial wave L. The phase shifts are directly related to the elastic scattering 
amplitude fL,d{p), where 
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In addition to having finite range, we assume also that the interaction is not too singular at 
short distances. Specifically, we require that the effective range expansion defined below in 



converges for sufficiently small p and that -^u'l\ is finite and u'l'^ vanishes as r — )• 0. 



,(p) 



Eq. 

For example, these short- distance regularity conditions are satisfied for a local potential, 
W{r^r') = V{r)d{r — r'), provided that V{r) = 0{r~'^~^'^) as r — )■ for positive e |30(]. In our 
discussion, however, we make no assumption that the interactions arise from a local potential. 
The treatment of spin-dependent interactions with partial wave mixing is beyond the scope 
of this analysis. For coupled-channel dynamics without partial wave mixing the analysis 
can proceed by first integrating out higher-energy contributions to produce a single- channel 
effective interaction. In order to satisfy our condition of energy-independent interactions, 
this should proceed using a technique such as the method of unitary transformation described 



uia pre 
in Ref . [silfi . 

The effective range expansion is 
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The term (5(dmod2),o is for odd d and 1 for even d. aL,d is the scattering parameter, ri^a 
is the effective range parameter, and are the n*'^-order shape parameters, pi^d is an 
arbitrary length scale that can be scaled to any nonzero value. The rescaling results in a 
shift of the dimensionless coefficient of p^^"*""^"^ on the right-hand of Eq. (|5]), and we define 
PL,d as the special value for pi^d where this coefficient is zero. 

Let u^l\ and w^j be radial solutions of the Schrodinger equation for two different mo- 
menta. We construct the Wronskian of the two solutions, 



(p) _^ (p') 

'^L,d^^^L,d 



ip') 
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and evaluate at some radius r > R. Taking the limits p' — )■ and then p — )■ 0, we find that 
for any r > R, 



21im 



dr' 



uld[r 
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where hi^dif) is defined as follows. For d = 2 and d = A, bi^dif) can contain logarithmic 
terms analog to Eq. ([5]). For the special case 2L + d = 2, we have 
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where 7 is the Euler-Mascheroni constant and for 2L + d = 4, 
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For the generic case of 2L + d any positive odd integer or any even integer > 6: 

2r(L + f - 2)r(L + f - 1) /r\ -2L-d+4 



IT 



1 



L + 



+ 



laL,d V2/ 

27T 



TiL + mL + l + l)al,\2J 



\ 2L~\-d 



(10) 



The formula in Eq. ( ITOj) for L = in three dimensions was first derived by Bethe [25j and 
extended by Madsen for general L (sij. The results presented here give the generalization 
to arbitrary d and arbitrary L. 

Since the integrand in Eq. ([7]) is positive semi-definite, ^ satisfies the upper bound 



rL,d < hAr) 



(11) 



for any r > R. For d = 3 our results are equivalent to the causality bound derived by 
Wigner j^. As noted in the introduction, the time delay between the incoming wave and 
the scattered outgoing wave is proportional to the energy derivative of the elastic phase 
shift, ddi^d/dE. Since the incoming wave must first enter the interaction region before 
the scattered wave can exit, causality places a bound on ddi^d/dE. The precise quantum 
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FIG. 1: Plot of 6L^3(r) — r^^s as a function of r for alpha-neutron scattering in the 5'i/2) -fi/2) ^-iid 
P3/2 channels. Causality requires this function to be non-negative for r > R. 



mechanical statement of this causality requirement is that the reciprocal logarithmic deriva- 
tive u"l\/ \ has a non-negative energy derivative. This fact can be derived from the 



Wronskian in Eq. ([6]), and a detailed derivation of this connection will be given in Ref. |35 



For finite-range interactions p2L+d-2 ^^^1- ^^^^ j^^^g ^ convergent effective range expansion, and 
d^L,d/dE at zero energy is proportional to the effective range r^ ^^. For d = 3, the Wigner 
causality bound in zero-energy limit is equivalent to the bound in Eq. ( ITTj) on the effective 
range. For S'-wave interactions in three dimensions the upper bound on the effective range 



was discussed in Ref. [26[. It was observed that for fixed ai^d the zero-range limit i? — )■ is 
possible only when tl^cL is negative. The constraint becomes more severe for larger 2L + d. 
For 2L + d > 4, the limit i? — )■ at fixed a^^d produces a divergence in the effective range, 
rL4 < hL,d{R) -00. 

Our results are exact only for the case where the interaction vanishes for r > R. For 
exponentially-bounded interactions of 0{e~^^^) at large distances, the results should still 
be accurate with only exponentially small corrections. For an exponentially-bounded but 
otherwise unknown interaction, the non-negativity condition for &L,d(r) — r^^d can be used 
to determine the minimum value for R consistent with causality. As an example, we plot 
^L,3(^) ~ fL,3 for alpha-neutron scattering in Fig. [T] In the plot, we show results for the 
'S'1/2, -Pi/2, and P3/2 channels. We note that a qualitatively similar plot was introduced 
for nucleon-nucleon scattering in the S'-wave spin-singlet channel j36|. The non-negativity 
condition gives R > 1.1 fm for S'1/2, R > 2.6 fm for P1/2, and R > 2.1 fm for P3/2. For 
comparison, the alpha root-mean-square radius and pion Compton wavelength are both 
about 1.5 fm. Since the minimum values for R are not small when compared with these, 
our analysis suggests some caution when choosing the cutoff scale for an effective theory of 
alpha-neutron interactions. 

At this point we comment on our requirement that the interactions are energy indepen- 
dent. For energy-dependent interactions it possible to generate any energy dependence for 
the elastic phase shifts even when the interaction W{r,r'; E) vanishes beyond some finite 
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radius R for all E. Under these more general conditions there are no longer any Wigner 
bounds and the constraints of causality seem to disappear. However, it is misleading to 
regard interactions of this more general type as having finite range. As noted in the intro- 
duction, the scattering time delay is given by the energy derivative of the phase shift. The 
energy dependence of the interaction can by itself generate large negative time delays and 
thereby reproduce the scattering of long-range interactions. In this sense the range of the 
interaction as observed in scattering is set by the dependence of W{r,r'; E) on the radial 
coordinates r, r' as well as the energy E. In this case the bound in Eq. ffTTl) can be viewed 
as an estimate for the minimum value of this interaction range. 

We now consider the scattering amplitude in the low-energy limit p — ?■ while keeping the 
interaction range R fixed. In the low-energy limit the scattering amplitude depends on just 
one dimensionful parameter when 2L + d < 3. For 2L + d = 1 and 2L + d = 3 the relevant 
parameter is ctL^d, and for 2L + ci = 2 it is pi^d- When 2L + > 4 a second dimensionful 
parameter appears in the non-perturbative low-energy limit. In the limit la^^dl ~^ oo, the 
upper bounds on the effective range reduce to the form p^ ii < ^e^ for 2L + d = 4, and 

r„ < _ 2r(L + i-2)r(L + f-l) ^^^^ 

for 2L + d > 5. There is no way to suppress the p'^\n{ppLd) term in the effective range 
expansion for 2L + d = 4 by fine-tuning parameters because the bound forbids tuning the 
argument of the logarithm to 1 as p — )■ 0. Similarly, for 2L + d > 5 the upper bound in 
Eq. f|T2|) and the negative coefficient on the right-hand side prevent setting ri^d to zero to 
eliminate the term ^r^^^^p^. Hence in both cases we are left with two relevant parameters in 
the non-perturbative low-energy limit. This corresponds to two relevant directions near a 
fixed point of the renormalization group, and the universal behavior is characterized by two 
low-energy parameters. For the case of P-wave neutron-alpha scattering in three dimensions, 
this issue was already discussed in [20]. Proper renormalization of an effective field theory 
for P-wave scattering requires the inclusion of field operators for the scattering volume 
and the effective range at leading order. In the renormalization group study of [l^, the 
emergence of multiple relevant directions around a fixed point was observed for the repulsive 
inverse square potential. 

For 2L + d > 4 this second dimensionful parameter has a simple physical interpretation 
for shallow bound states. Consider a bound state ai p = ipi in the zero binding-energy limit 
P/ — )• 0^. Let P>(r) be the probability of finding the constituent particles with separation 
larger than r: 



P^(r) = I dr' 



oo 



(13) 



where u^l^j"^ is the normalized wave function. For 2L + d < 3, the probability P>(r) equals 1 
in the limit pi — ?• O"*" for any r. At sufficiently low energies the physics at short distances is 
irrelevant, and the bound state wavefunction is spread over large distances. For 2L + d > 4, 
however, the situation is different. For 2L + d = 4 the probability is logarithmically 
dependent on p^ ^ and can be tuned to any value between and 1 |35|. Similarly for 
2L + d>5, 

2r(L + | 2)r(L + i-i) 

for r > R. For this case the characteristic size of the bound state wavefunction is 
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For P-wave Feshbach resonances in alkali atoms our analysis must be modified to take 
into account long-range van der Waals interactions of the type W{r,r') = —CQr~^6{r — r') 
for r, r' > R. This raises various new issues such as the applicability of our approach 
to power law potentials and the appearance of non-analytic terms in the effective range 
expansion. These issues will be addressed in detail in Here we will briefly discuss 



the modifications for potentials with a van der Waals tail in three dimensions only. It is 
convenient to reexpress Cq in terms of the length scale (3^ = {2ixCeY^'^. In the following, we 
set (i = 3 and drop the d subscript. Instead of free Bessel functions, scattering states should 



be compared with exact solutions of the attractive r ^ potential [37|, |38j. The effect of the 



interactions for r < R are described by a finite-range i^-matrix Kl{p^) which is analytic 
in [39I, Kiip^) = ^n=o 1 - -^l"^ P"^" ■ When phase shifts are measured relative to free 
spherical Bessel functions, the effective range expansion is no longer analytic in p^. For 
L = 0, the leading non-analytic term is proportional to p"^. For L = 1 the non-analytic term 
is proportional to p^, thereby voiding the usual definition of the effective range parameter. 

For a pure van der Waals tail, however, one can still obtain useful information from our 
approach. The zero-energy resonance limit is reached by tuning the lowest-order A'-matrix 
coefficient K^^ to zero. It turns out that for L = 1 in this limit the p^ coefficient in 
the effective range expansion also vanishes, and we can therefore define an effective range 



parameter for both S- and P-waves [38|, |40 



ro = [r(l/4)f (/36 + 3irf/3e-i) /(3vr), 

n = -36 [F (3/4)]^ - 5K?^^^') /(Stt). (15) 

(2) 

For the case of single-channel scattering for alkali atoms, the coefficients are negligible 

compared with This is also true for some multi-channel Feshbach resonance systems 

(4l| . In these cases we observe that the upper bounds for r^, in Eq. (IT^ are satisfied for 

L = and L = 1 when we naively take P ~ /Sg- In general, there may be multi-channel 

(2) (2) 
systems where the coefficients cannot be neglected. Nevertheless, the coefficients K^^ 

should satisfy Wigner bounds similar to those derived here for the effective range. This may 

be a useful starting point for further investigations of multi-channel Feshbach resonances in 

alkali atoms. 

In this letter, we have addressed the question of universality and the constraints of causal- 
ity for arbitrary dimension d and arbitrary angular momentum L. For finite-range inter- 
actions we have shown that causal wave propagation can have significant consequences for 
low-energy universality and scale invariance. We find that in certain cases two relevant low- 
energy parameters are required in the non-perturbative low-energy limit. In the language 
of the renormalization group, this corresponds to two relevant directions in the vicinity of 
a fixed point. In particular, we confirm earlier findings in the case three dimensions for 
P-wave scattering 2^ based on renormalization arguments and for higher partial waves in 



general [2J] in the framework of the renormalization group. The analysis presented here 
concerns only the question of universality in two-body scattering. Universality for higher 
few-body systems requires a detailed analysis for each system under consideration. Effective 
field theory and renormalization group methods may again provide a useful starting point 
here Our results may help to clarify some of the conceptual and calculational is- 



sues relevant to few-body systems for general dimension and angular momentum and their 
simulation using short-range interactions. 
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